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Abstract 

Using 24.5 million ip(2S) decays collected with the CLEO-c detector at CESR we present the 
most precise measurements of magnetic dipole transitions in the charmonium system. We measure 
B(i/j{2S) -» 7%) = (4.32±0.16±0.60) x 1(T 3 , B(J/ip -» ^r] c )/B{tp(2S) -» 7%) = 4.59±0.23±0.64, 
and B{J/ij} — > 7%) = (1.98 ±0.09 ±0.30)%. We observe a distortion in the r/ c line shape due to the 
photon-energy dependence of the magnetic dipole transition rate. We find that measurements of 
the 77c mass are sensitive to the line shape, suggesting an explanation for the discrepancy between 
measurements of the rj c mass in radiative transitions and other production mechanisms. 
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The spectrum of bound charm quarks provides an important testing ground for our un- 
derstanding of Quantum Chromodynamics (QCD) in the relativistic and non-perturbative 
regimes. Radiative transitions in particular have recently been the subject of both lattice 
QCD calculations [lj and effective field theory techniques [2]. Key among these are the 
magnetic dipole (Ml) transitions J/tp — > 7% and ip{2S) — > 77/c which are among the most 
poorly measured transitions in the charmonium system. Not only are precision measure- 
ments needed to validate our theoretical understanding, precise measurements of these Ml 
transitions are critical for normalizing rj c branching fractions, a key input to extracting other 
properties such as r 77 (^ c ). The J/ip and ip(2S) — > •yr] c transitions are also a source of in- 
formation on the 7] c mass and width. There is currently a 3.3 a inconsistency in previous r) c 
mass measurements from J/ip and ip(2S) — > jrjc (averaging 2977.3 ± 1.3MeV/c 2 ) compared 
to 77 or pp production (averaging 2982.6 ± 1.0MeV/c 2 ) [3j. 

In this Letter, we present the most precise measurements of B(J/ip — > 7^) (abbrevi- 
ated Bis), S('0(25') — > 7^ c ) (abbreviated Bos), an d their ratio using 24.5 million ip(2S) 
decays collected with the CLEO-c detector [4j. For the first time, we clearly observe the 
distortion of the rj c line shape in the photon energy spectrum due to phase space and energy- 
dependent terms in the Ml transition matrix element. We find that this line shape distortion 
may be responsible for the inconsistency in measured r] c mass. 

The CLEO-c detector operates at the Cornell Electron Storage Ring (CESR) which 
provided symmetric e + e~ collisions at the ip{2S) center-of-mass. The detector features a 
solid angle coverage of 93% for charged and neutral particles. The charged particle tracking 
system operates in a 1.0 T magnetic field along the beam axis and achieves a momentum 
resolution of ~0.6% at p = 1 GeV/c. The cesium iodide (Csl) calorimeter attains photon 
energy resolutions of 2.2% at = 1 GeV and 5% at 100 MeV. Two particle identification 
systems, one based on ionization energy loss (dE/dx) in the drift chamber and the other 
a ring-imaging Cherenkov (RICH) detector, are used together to separate from tt^. 
Detection efficiencies are determined using a GEANT-based [fl Monte Carlo (MC) detector 
simulation. 

To extract B2S we use the ~ 640 MeV photon transition line visible in the inclusive photon 
energy spectrum from multi-hadronic events collected at the ip{2S) resonance. A series of 
exclusive decay modes of the rj c (where the background can be greatly suppressed) are used 
to constrain the line shape for the inclusive spectrum. To measure Bis/Bis, we take the 
ratio of events in the chains 

ij(2S) -> 7T+7T- J/fr J/iJ> -> 1Vc ; rj c -> X t (1) 

V(25) - ^Vc, Vc - X h (2) 

where the Xi are exclusive decay modes of the i] c ; we then adjust for efficiencies and 
B(ip{2S) — > 7r + 7r~ J /ip). Rather than using the inclusive photon spectrum from J /if) de- 
cays, we minimize systematic errors by taking B\s to be the product of Bis with B\sjB<is- 

The first half of this Letter describes three samples of events: exclusive decays of the J / ip 
(Eq. flTJ); exclusive decays of the ip(2S) (Eq. ([2])); and the inclusive photon spectrum from 
ip{2S) decays. We use the exclusive samples to investigate the photon-energy dependence of 
the rj c line shape. In the second half, our measurement techniques, guided by our line shape 
investigations, are more fully developed. 

Twelve exclusive rj c decay modes are used: X^ = 2(7r + 7r~), 3(7r + 7r~), 2(7r + 7r~7r°), 
nT^Ks, tc°K + K~, 7i^7i + 7i- ^Kg, n + 7T-7r°K + K~, 2{K + K~), 2{n + n~)K + K~ , 
7i + 7r~K + K~ , 7r + 7r~?7, and 2(n + 7r~)ri, where the r\ is detected in either its 77 or 7r + 7r~7r° 
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decay mode. These include all previously reported decay modes of the r\ c (except pp, which 
has a comparatively small rate) 0] in addition to new decay modes observed here with 
comparable raw yields. 

The reconstruction of the J /if} and if)(2S) — > 7% exclusive decay chains share several 
selection criteria. In addition to standard fiducial requirements, photons must have energy 
greater than 30MeV and must not align with the projection of any track into the calorimeter. 
For 7T° and rj decays to 77, the mass of the pair of daughter photons is required to be within 
3 a of the nominal mass. We require fitted tracks of charged particles to have xVd.o.f. < 50 
and be located in the central region of the detector (| cos#| < 0.93, where 9 is the polar 
angle with respect to the e + direction). All charged tracks must be positively identified by a 
combination of dE/dx and the RICH detector. To reconstruct 77 — ► 7r + 7r~7r° all three decay 
products must pass the above criteria and must have an invariant mass within 30MeV/c 2 
of the nominal 77 mass. The Kg candidates are selected from pairs of oppositely charged 
and vertex-constrained tracks with invariant mass within 15 MeV/c 2 of the K% mass. A 
four-constraint kinematic fit of all identified particles to the initial ip(2S) four-momentum is 
performed and a x| c /d.o.f < 5 is required. This both sharpens the measured momenta and 
reduces backgrounds due to missing particles or particle misidentification. The hypothesis 
with the best fit quality is accepted per decay mode; less than 0.5% of these events enter 
multiple modes. 

For the selection of exclusive J /if) decays, the recoil mass of the 7r + 7r~ pair is used to 
select the process if>(2S) — > n + n^J/ip and to select J /if) sidebands. A sideband subtraction 
is used to account for non-J/ip decays. The 7r + 7r~ recoil momentum is used to boost the 
photon energy into the J /if) rest frame. 

Fits to the resulting photon energy spectrum for the sum of all r\ c decay modes are shown 
in Fig. [H The background shape has two essential features: (i) background that falls with 
energy from J /if) — > Xi, where a spurious cluster is found in the calorimeter, and whose 
shape is modeled by MC simulation [3] and (ii) a rising background from both J /if) — > TT°Xi 
and non-signal J /if) — > 7 A; that is freely fit to a second degree polynomial. (The polynomial 
form is motivated by MC and is validated by comparing reconstructed J /if) — > 7r° Aj in both 
data and MC.) These two contributions and their total are shown as dot-dashed and dashed 
curves. A fit using an unmodified relativistic Breit-Wigner distribution (dotted line), with 
the amplitude, mass, and width as free parameters, fails on both the low and high sides of 
the signal. A fit using a relativistic Breit-Wigner distribution modified by a factor of j2| 
improves the fit around the peak but leads to a diverging tail at higher energies (not shown). 
To damp the E^ an additional factor of exp(— E^/f3 2 ) is added, inspired by the overlap of 
two ground state wave functions. The resulting fit has a 25% confidence level (solid line), 
with f3 = 65.0 ± 2.5 MeV. In all cases the signal shapes used are convolved with a resolution 
function determined from MC simulation. The resolution is 4.8 MeV after the kinematic fit. 

The uncertainty associated with the line shape prohibits precision mass and width mea- 
surements. It is interesting to note, however, that the resulting t] c mass in the unmodified 
Breit-Wigner fit is 2976.7 ± 0.6 MeV/c 2 (statistical error only), consistent with previous 
measurements from J /if) and ip{2S) — > 777 c , while the modified Breit-Wigner fit returns an 
T) c mass of 2982.2 ± 0.6 MeV/c 2 (statistical error only), consistent with that determined from 
77 fusion and pp annihilation. To resolve this inconsistency, a thorough understanding of 
the r] c line shape in J/ if) and if)(2S) — ► 77/ c will be required. 

For the selection of exclusive ip(2S) decays, the transition if>(2S) — > n+n^J/ip is sup- 
pressed by requiring there be no pair of oppositely charged particles (assumed to be pi- 
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ons) with recoil mass within 15 MeV/c 2 of the J/i/> mass. The photon energy spectra 
for individual r] c decay modes are shown in Fig. [2j the sum of all modes is shown in 
Fig. [3K Several small nonlinear backgrounds below 560 MeV are apparent and are due 
to a combination of (i) if)(2S) — > n°h c ;h c — > jt] c ] (ii) ip{2S) — > '"fXcj', Xcj jJ/i^] an d 
(hi) ip(2S) — > ir°J/il). Based on detailed MC studies, all other backgrounds are linear, the 
largest being ip{2S) — > 7r°Xj. 

Fits to the ip{2S) — » 777,- photon energy spectrum with a relativistic Breit-Wigner dis- 
tribution convolved with an experimental resolution function (with a resolution of 5.1 MeV 
after the kinematic fit) were unsuccessful. For a hindered Ml transition the matrix element 
acquires terms proportional to E^, which, when combined with the usual E^ term for the 
allowed transitions, lead to contributions in the radiative width proportional to E'L [2]. We 
find that if we assume a linear background, as indicated by MC simulations, we are not able 
to obtain a good fit to our E 1 spectrum for the sum of exclusive ip(2S) — > 7% modes with a 
pure E^ dependence. We therefore use the empirical procedure described below to extract 
the ip(2S) — > jr/c yield. 

Extensive cross-checks have been performed to verify that the line shape asymmetry is 
not an experimental artifact. Events selected without the aid of a kinematic fit indicate 
an asymmetric line shape independently in both the photon energy and the hadronic mass. 
The asymmetric line shape is not correlated with r] c decay modes that include 7r°, Kg, or 
77 candidates. No indication of either asymmetry or peaking background has been found in 
detailed MC studies, where all known decays in the charmonium and light quark systems 
are simulated and unknown decays are modeled with the EvtGen generator Q]. The photon 
angular distribution from ip(2S) —> jrj c fits the 1 + cos 2 9 distribution expected for Ml tran- 
sitions, whether using symmetric or asymmetric signal shapes to extract yields in different 
regions of cos 9. 

For the selection of the third sample of events, the inclusive photon spectrum from ip(2S) 
decays, the photon is required to pass the same requirements and ip(2S) — ► 7T + 7T~J/ifj 
is suppressed in the same manner as the exclusive ip(2S) — ► 7?7 C . To suppress the 7r° 
background, each signal photon candidate is paired with all other photons in the event and 
is rejected if the invariant mass of the pair is within three standard deviations of the 7r° 
mass. Backgrounds due to e + e~ QED processes are substantially reduced by requiring one 
or more charged tracks in an event in combination with requirements on total energy and 
track momenta j^]. The inclusive photon spectrum is shown in Fig. [3fc. The rise at lower 
energies is from Xcj jJ/i*- All other backgrounds are smooth and are dominated by 7r° 
decays. 

Our final results are obtained from: 



N INC 

^25 = inCtJ " ( 3 ) 



Bis Nf s xc 



(4) 



" 1S = _INC /_EXC /_EXC\ at K! 



B2S N§<c (^c/efxc) ^ 
ego (eW C /eW C )^(2S)B m 



where N and e represent the observed yields and the calculated efficiencies of ip(2S) and 
J/ijj in inclusive (INC) and exclusive (EXC) rj c channels. The B(ip(2S) — > 7t + tt~ J/ijj), 
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TABLE I: Final yields and efficiencies. 



N^ c 59510 ± 2145 

Nf s xc 5376 ± 199 

N 2S C / N fs XC 11-07 ±0.33 
Nf s xc 5638 ± 187 

4^ c 56.37% 
£iS C /4s C 0-6515 



TABLE II: Summary of systematic errors (in percent). 



Systematic Error (%) B2S Bis Bis/&2S 

Fitting for N^ c 8 - - 

Fitting for Nf,? 8 

Fitting for N^ c /Nf s xc 4 

Fitting for Nff 10 10 

Effect of Line Shape on 4s° and ef s xc 7 4 3 

Effect of Line Shape on xc 1 1 
MC Modeling of Inclusive rj c Decays 

and Inclusive Event Selection 8 8 

Exclusive Efficiency Ratio Calculation 2 2 

Exclusive Event Selection - 3 3 

640 MeV Photon Efficiency 2 2 

110 MeV Photon Efficiency 2 2 

tt + it~ Efficiency - 2 2 

^ip(2S) 2 2 

B(ip(2S) TT+vr- J/VO 2 2 

Total 14 15 14 



abbreviated B nn , is taken from a previous CLEO measurement, (35.04 ± 0.07± 0.77)% [loj ]. 
N^/,(2S) is the number of ^p(2S) decays, 24.5 million, which is known to 2% [lOj]. Final values 
are listed in Table [I] Systematic errors are listed in Table HT1 

The procedure used to obtain N™ , Nif , and N^ c /Nf^ c is threefold: (i) the back- 
ground to the exclusive ip{2S) — > jt] c process is fit with a first order polynomial using 
regions above and below the signal (Fig. [3^); (ii) the background is then fixed and carried to 
the exclusive if>(2S) — > ^rjc spectrum that has not been adjusted by a kinematic fit (Fig. [3b), 
which is directly comparable to the inclusive spectrum; (iii) the histogram obtained by sub- 
tracting off the background in (ii) is used to fit (along with a fourth degree polynomial for 
the background) the inclusive photon spectrum (Fig. [3fc). Numbers of events above back- 
ground are counted for photon energies between 560 MeV and 1100 MeV. As can be seen 
from the background subtracted inclusive spectrum (Fig. [3]i) , there is excellent agreement 
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between the exclusively determined signal shape (line) and the signal shape present in the 
inclusive photon spectrum. 

Systematic errors due to this fit procedure are determined by varying the orders of the 
background polynomials, the ranges of the fits, and the range of the signal that is excluded 
in the exclusive fit. We find 8% variations for N™ and Nf^ , but only 4% variations in 
their ratio. In addition, since we integrate our signal between 560 MeV and 1100 MeV, our 
uncertainty in the line shape affects the signal efficiency We take a conservative systematic 
error of 7% for e J ^ c to cover the two extreme cases that the signal shape has a Breit-Wigner 
tail and the case that the entire signal lands within the signal region. The systematic error 
is smaller for efg (3%) because the kinematic fit pulls more of the signal within the signal 
region. These errors are correlated resulting in an error of 4% in the ratio. 

The measurement of N^ is from the modified Breit-Wigner fit shown in Fig. (TJ which 
results in an r] c width of 31.5 ± 1.5 MeV/ c 2 (statistical error only). A fit using an unmodified 
Breit-Wigner distribution with the i] c width fixed to the current PDG value (26.5 MeV/c 2 ) 
gives a smaller value of Nf^ by 10%. A 10% systematic error covers this extreme minimum, 
as well as variations using a Breit-Wigner distribution modified only by an term and 
variations of the background parametrization. The signal is integrated above 40 MeV, which 
introduces an additional 1% systematic error into the efficiency estimate (efg^ ) due to signal 
shape uncertainty. All efficiences are independent of photon energy in the signal regions; 
hence the signal shape uncertainty only affects efficiencies by way of the limits of integration. 

Because only ^25% of all r\ c decays are known, our estimate of depends on modeling 
the unknown decays of the r] c . We use two different hadronization models, but then weight 
the efficiencies for different track multiplicities according to multiplicities observed in data. 
We account for our limited sensitivity to events with no tracks by varying their weight by 
±100%. We also simultaneously vary the event selection requirements, keeping and removing 
the 7T° suppression, and using several QED suppression schemes. All variations are covered 
with an 8% systematic error. 

The ratio of exclusive efficiencies, sfs / £ 2S' * s calculated by weighting the efficiency 
ratio for each rj c decay mode individually using the number of ip{2S) — > / ~fi] c events observed 
in each mode. The ratio of efficiencies has a slight dependence on decay mode. Varying the 
number of events in each mode by 30%, we find a 2% systematic error due to the composition 
of decay modes. 

Many systematic errors, such as those due to tracking or 7r° reconstruction efficiencies, 
cancel in the ratio of exclusive efficiencies, £fs/ £ fs ■ To estimate any possible dependence 
of the final numbers on the exclusive event selection, a wide range of requirements on the 
quality of the kinematic fit are imposed, and photons are required to be in different regions 
of the calorimeter, resulting in variations of less than 3%. We conservatively assume that 
systematic uncertainties for the reconstruction efficiencies of the ~ 110 MeV and ~ 640 MeV 
transition photons do not cancel in the ratio and assign 2% uncertainties for each. We also 
assign a 2% uncertainty for the efficiency of the 7t + tt~ from ip(2S) — > -k + -k~ J/ip. 

We find B 2S = (4.32±0.16±0.60) x 10~ 3 , B 1S /B 2S = 4.59±0.23±0.64, and B 1S = (1.98± 
0.09±0.30)%. Both Ml transitions Bis and B 2 s are larger than previous measurements (9l.Hl| 
(56% and 44% higher than the PDG averages [3j, respectively) due to a combination of a 
larger rj c width and an accounting for the asymmetry in the line shapes. 

In conclusion, we have studied J/ip and ip(2S) — > 7// c transitions and measured their 
branching fractions. These new measurements will renormalize many rj c branching fractions 
and T 77 (?7 C ). We also find that a thorough theoretical understanding of the r\ c line shape in 
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Ml transitions in the charmonium system will be crucial if the systematic errors on these 
branching fractions are to be reduced and if the mass and width of the r\ c are to be extracted 
from these processes. 
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FIG. 1: Fits to the photon spectrum in exclusive J/tp — > jrj c decays using relativistic Breit- 
Wigner (dotted) and modified (solid) signal line shapes convolved with a 4.8 MeV wide resolution 
function. Total background is given by the dashed line. The dot-dashed curves indicate two major 
background components described in the text. 
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FIG. 2: The kinematically fitted photon energy spectra from ip(2S) — > ^J] c for individual r/ c decay 
modes. 
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FIG. 3: (a) The kinematically fitted photon spectrum from the sum of exclusive ip{2S) — > "fq c 
modes with a polynomial fit to the background (solid line for the regions included in the fit; dotted 
line elsewhere), (b) The photon energy spectrum that has not been adjusted by the kinematic fit 
with the same background shape as (a) overlaid, (c) The fit to the inclusive photon spectrum in 
ip(2S) decay. The signal shape (solid line) is described in the text. The background is given by 
the dashed line, (d) The background subtracted inclusive photon spectrum with the signal shape 
overlaid (solid line). 
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